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Abstract 

We consider travelling periodical and quasiperiodical waves in sin- 
gle mode fibres, with weak birefringence and under the action of cross- 
phase modulation. The problem is reduced to the "1:2:1" integrable 
case of the two-particle quartic potential. A general approach for find- 
ing elliptic solutions is given. New solutions which are associated with 
two-gap Treibich-Verdier potentials are found. General quasiperiodic 
solutions are given in terms of two dimensional theta functions with 
explicit expressions for frequencies in terms of theta constants. The 
reduction of quasiperiodic solutions to elliptic functions is discussed. 
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1 Introduction 



We consider the system of two coupled nonlinear Schrodinger equations 



iAz + -Att ^ oB ^ {\Af + -i\B\^)A = 0, 
iBz + ^Btt + aA+ + j\A\^)B = 0. 



(1.1) 



These equations are important for a number of physical applications. For ex- 
ample the dimensionless circularly polarized components A, B of the electric 
field in a single-mode straight (no twists) fibre satisfy these equations ^ 



If the functions A and B do not depend on the variable t, these equations are 
the generalized discrete self-trapping dimer system [|l^, which are integrable 
in terms of elliptic functions. If 7 = these equations are the well known 
DST dimer equations [1^. The equations (|1 . 1|) do not belong to the class 



of nonlinear evolution equations integrable by means of the inverse scatter- 
ing method |2^. Nevertheless, they can have exact vector solitons P, 3C], 
bound solitary waves M and periodical solutions |T§, EDI. This interesting 



phenomenon is related to the existence of a fourth integral found by Dowling 
Due to this property, two-soliton solutions also exist ^0 . 



In the present paper the general polarization-modulated states - two- 
gap solutions - are expressed in terms of two-dimensional theta functions. 
Another motivation for this work is the attempt to classify the exact solutions 
of the Coupled Nonlinear Schrodinger equation according to the program of 
classification given in 

Let us introduce the new functions A = {a + ib)/ ^/2, B = {a — ib)/ y/2 
and new independent variables z = + 1) Z/2, t = y/'j + IT |jl6|. Then the 
equations (|1 . 1| ) can be rewritten as follows 



ia^ + au + QqG + p{\a\'^ + \b\'^)a + q{a^ + b'^)a* = 0, 
lb, + btt - nob + p{\b\^ + \a\^)b + q{a^ + b^)b* = 0, 

where Qq = 2cr/3, p = cr/3 and q = 1/3. 

We seek solutions of equations (|1.2|) in the following form W6 



[1.2) 



a{z,t) = qi{t) exp{iQz) , b{z,t) = q2{t) exp{iQz), 



[1.3) 
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where qi{t) and q2{t) are real functions, and Q is an arbitrary real constant. 
We obtain the system 



(^2 

+ {ql + ql)q2 = + ^o)q2. (1-4) 

These equations describe the known integrable case "1:2:1" of the quartic 
potential W = Agf + _Bg^g| + Cg| [[l^, which is separable in ellipsoidal coordi- 



nates. The foregoing analysis can be applied also to the integrable potential 
"1:12:16" which is separable in parabolic coordinates. The dynamics of the 
other separable cases of the quartic potentials "1:6:1" and "1:6:8" can be 
expressed in terms of two elliptic functions with different moduli [p7 |. 



The paper is organised as follows. In Section 2 we describe the Poisson 



structure of the integrable system "1:2:1" using the results given in ITT], [12 
We derive a Lax representation for the system, prove its complete integrabil- 
ity and construct the separated coordinates. In Section 3 we show that the 
problem of the description of elliptic solutions for the system ( |1.4] ) is reduced 
to the description of the elliptic potentials for the Schrodinger equation and 
the construction of the corresponding Lame polynomials. The problem is 
closely related to the structure of the locus for the Calogero-Moser system 
[0]. In Section 4 we analyse the known elliptic solutions for the system (|1.4| ) 
and show that they are included into the approach developed. Moreover we 
find a new elliptic solution for the system ( [L.4D which is associated with the 
two gap Treibich-Verdier potential |^1|. In Section 5 we give an integration of 



the system via theta functions and calculate the frequencies in terms of theta 
constants. The question of the periodic solutions is formulated from a theta 
functional point of view and reduced to the description of some modular va- 
rieties given by the vanishing of some two dimensional theta constants. Such 
a reduction is demonstrated on the example of the elliptic solution associated 
with the Treibich-Verdier potential. 
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2 The Lax representation 



The system (L4) is a completely integrable Hamiltonian system with the 
Hamiltonian 

H = lipl+pl) + \iql + qlf - \{z,ql + z,ql), (2.1) 

where {qi,Pi), 2 = 1,2 are canonical variables with Pi = (ji = dqi/dx, i = 1,2, 
and Zi,Z2 are arbitrary constants. 

The Lax representation for the system ( p.. 41) is known in terms of 3 x 3 
matrices |Q . However, we find it convenient to use a 2 x 2 Lax representation, 
L{z) = [L{z), M{z)] which is a special case of the hierarchy of separable 



systems discussed in [0, [12 



To do this we fix the following ansatz for the Lax operator. 



L{z) 

where 



^ M(z) = ( ^ ^] (2 2) 



U{z) = Viz) = ~lu{z), (2.3) 

= ^E-37-^(^)' Q^^) = ^AjZql (2.4) 

The associated algebraic curve det(L(z) — yl) = has genus two, and is 
written as 

= -4{z-zi){z-Z2){z'^ -z\zi + Z2)+z{ziZ2-H)-F), (2.5) 
= -Alliz-z,), 

i=l 

where w = y{z — zi){z — Z2), H is the Hamiltonian (^.l]), F is the additional 
integral of motion, 

F = -^{Piq2 - P2qi) + 2 (?i + ^2) [Z1Z2 - —qi - — ^2 

-\{z2pI-z,pI), (2.6) 
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and 2:3, Z4,, are the roots of the cubic on the rhs. of (p.5|). Let us define new 
coordinates /ii,/i2 as zeros of the function U{z) in the Lax operator, i.e. 
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- /Ui)(z2 - /U2) 



(2.7) 



On can prove (see e.g. 
are defined as 



Zi — Z2 Z2 — Z\ 

TT| , [T2[) that the canonically conjugated momenta 



Wi 



vr2 



Z2 



/i2 - /il 



(2.8) 
(2.9) 



2 (/ii - zi)(/ii - 2:2) 2 (/i2 - 2:1) (Ai2 - Z2) 

and therefore the coordinates (/ii, tTj) are the separated coordinates (eUip- 
soidal coordinates). 

It follows from (|2.7| , |2.9|) that the the dynamics of the system described in 
the coordinates [ni, fii) become the Jacobi inversion problem associated with 
the curve (12. f 



Mo 



W{fi) 



+ 



M2 



MO 



w(/i) 



Mo 



w(/i) 



+ 



Mo 



w(/i) 



2x + 6, 



(2.10) 



where a, 6 are constants defined by the initial conditions. If zi^2 are real and 
zi > Z2 then the solution ([2^^ ) is real if /i2 < -^i < yUi and yUi < Z2 or Z2 < fi2- 
We remark that the integrable case "1:12:16" of the quartic potential 
arises in this approach by the introduction of another ansatz for the function 
U whose zeros define parabolic coordinates W^. 



3 Periodic solutions associated with the Lame 
equation 



In this section we develop a method (see also JTO], |T^, |T9[) which allows us to 
construct periodic solutions of (|1.4| ) in a straightforward way. The method 
is based on the application of spectral theory for the Lame equation with 
elliptic potentials ||l|, 

— ^^(x, z) - W(x)^(x, z) = -z'ifix, z). (3.1) 
dx^ 
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with U (x) being an elliptic potential. 

Because the algebraic curve (|2.5|) associated with the problem has genus 
two we shall consider a two-gap elliptic potential for the equation Such 
potentials are known to be of the form [|^ 

N 

U{x)=2Y^p{x-Xi), (3.2) 

i=l 

where p{x) is the Weierstrass elliptic functions with the real period 2a; = 
2ljJi and imaginary period 2uj' = 2003. The number is a positive integer 
N > 2 (the number of "particles") and the numbers x = (/2i, . . . , xn) belongs 
to the locus Cn, i-e., the geometrical position of the points given by the 
equations 

£^ = I (x); 5] p'{xi - X,) = 0, J = 1, . . . iV)| . (3.3) 

Equation ( p.l|) allows the coalescence of three particles Xi and the potential 
takes the form 

n m 

W(x) =6^p(x-a;i) + 2^p(x-Xj), 3n + m = N (3.4) 

i=i j=i 

The associated algebraic curve of genus two can be described with the 
help of the Novikov equation For example, let us consider the two-gap 
potential for (|3.4| ) normalized by its expansion near the singular point as 
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U{x) = — + ax^ + bx^ + cx^ + dx^ + 0{x^^). (3.5) 



Then, following from the Novikov equation ||2^, the algebraic curve associ- 
ated with this potential has the form 0] 

= h \ UH . (3.6) 

2 2 \ 8 A j 4 2 ^ ' 

Let us consider the trace formulae written for the elliptic potential 
in the form 

N 5 

/^i+/i2 = -IZp(a; -a;j) + 

i=i ^ j=i 
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/il/i2 = 3^p{x - Xi)p{x - Xj) 



Ng2 



i<j 

+^E-^-.-^(e%) (3-7) 

We require that the eigenvalues Zi and Z2 be the branching points of the 
curve associated with the two-gap elhptic potential Using the relation which 
follows from ( |2.7| ) 

ql + ql = 2{fii + fi2 - zi - Z2), (3.8) 
then from the first trace formula we can write ( |1.4| ) in the form 

— gi - U{x)qi = {n-no + 2^1 + 2z2)qi (3.9) 

—q2 - U{x)q2 = (fi + ^]o + 2^1 + 2^2)52, (3.10) 

where we set without loss of generality Z]j=i = 0. 
Then if the relations 

Vt - VLq = -3zi - 2z2, Vt + VLq = -2zi - ?,Z2. (3.11) 

hold, the problem of finding elliptic solutions for ( |1.4| ) is reduced to the 
calculation of two Lame polynomials, i.e. the values of the eigenfunctions of 
( |2.2D corresponding to the spectral parameter fixed at the ends of the gaps 

gi = Ci^(x, ^i), q2 = C2^(x, Z2), (3.12) 

where Ci are some constants, corresponding to a point x on the locus Cjsi. 

For the elliptic potential without degeneracy ( p.2|) we have from ( |2.7|J3.71 ) 
and ( p.8|) the general formula for elliptic solutions of equations (|1.4|) . 



1 / ^ 
ql = 2zl + 22:1 J2 Pi^ - ^i) 



Ng2 



+ 6 ^ p{x - Xi)p{x - Xj) — + ^ ZiZj\ , (3.13) 

l<i<j<N l<i<j<5 



1 ( ^ 

2zl + 2z2 V p{x - Xi) 

Zl-Z2\ ^ 



Ng2 



+ 6 ^(^ - Xi)p{x - Xj) — + I • (3-14) 

l<i<j<Af ^ l<i<j<5 
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A problem with the apphcation of this formulae is that we have to find 
points X on the locus such that the functions ( p.lS , ^?!! ) are real and finite. 
Some special cases are known, for example the analytical description of the 
locus ( |3.3| ) for two-gap elliptic potentials for which Xi are half-periods, Xj = Ui 
(Treibich-Verdier potentials, 0]). Details are given in [|l^. The curves and 
Lame polynomial associated with these potentials are given in Table |I|. 



4 Elliptic solitary waves 

We shall show in this Section that the known periodic solutions of (|1.4| ) 
are associated with one and two-gap elliptic potentials of the Schrodinger 
equation. 

4.1 One-gap potentials 

The equations (|1.4| ) have the following solutions [^: 

qi = Cisn(ax, k), q2 = C2cn{ax, k), (4.1) 

where the amplitudes Ci ,6*2, the modulus k of the elliptic functions, and the 
temporal pulsewidth 1/a of the waves are defined in terms of the physical 
parameters Q and Qq as 

a^k^ = 2fio, Cf = n- 3^0 + a^ C| = ^] + fio + a^- 

The parameters can be expressed in terms of the Weierstrass parameters 63 < 
62 < ei as follows, = -(e2 + 63)72, ilo = (62 - 63)72, a = y/ei - 63, Ci = 
Vei - 262, C2 = Vei - 263 

The solutions ( [4.1| ) are associated with the eigenvalues zi = 62 and Z2 = 63 
of the one-gap Lame potential. In this case H = |(62 + 63)^ -|- 6262, F = 
1(62 + 63)^ and the curve ( p^) reduces (after a transformation of the spectral 
parameter z ^ — 2; + 62 + 63) to the product of the Weierstrass cubic and a 
perfect square 

/ 3 \2 

= 4{z - ei){z - e2){z - es) [z - -eij . (4.2) 

It is easy to see that one of the variables fii is pinched at the point 36i72 
and the Jacobi inversion problem written for the curve (^4.2| ) becomes the 
inversion of the elliptic integral for the remaining variable /i. 
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Table 1: The algebraic curves and Lame polynomials 





Tlie c-ur\'es and Lame poh'noniials A(,/-) 




U(x) — 6p(x) 

^ij ^ \/ipi^) - ei)(p(x) - Cj), (z^Sck), i 7^ /c = 1,2,3 

A± = p{x)±l^, (z = ±V3^) 




U{x) = 6p{x) + 2p{x + Ui) 
w^ = -{z + 6e,) Ut=iiz - Zk{i)), i = 1, 2, 3 
^1,2(0 = e,- + 2ei ± 2^ (e^ - ej)(2ej + 7ej) 
2^3,4(0 = Cfc + 2ei ± 2^/ (ci - efe)(2efc + Tcj) 


Aik^^{p{x) ei){p{x) ej) + l[{ei e^) ± ^J {ci ek){7ei + 2ek)]^ l\J^J^ 

{z ^ ek + 2ei ± 2^{ei - ek){2ek + 7ei)) 
Ao = p{x) -ei{z^ -Gcj) 


W5 


U{x) = Qp{x) + 2p{x + Ui) + 2p{x + tUj) 

w^ = mliz-z^ik)), J = 1,2,3 
Z4{k) = 6efc - 3ei, z^ik) = 6ei - 3efe 

nf=i(;^ - = - Sc'efc + ihlel - 2Qg2)z - 369e^ + 1326^^2 
K-sJ{p{x) ei){p{x) e,) + (e, e^O^fg^J, - 6e, 3e,-) 

ej)(p(x) ek) + {ej ei)^^J^J^_J:, (z - 6ej 36^) 


A„ = - e,)(p(x) - e,) + a.^yggf^ + ^..^gg^, 

f \ ■ _L ■ _L 1 1 Q ~ 15e?+27e|-6eiej-4+22„(ej-ei) 
1,-2 — Zn)-i ^ 7^ J 7= /C, n — i, Z, 0, — 24(ej-ei) 



10 



4.2 Periodic solutions associated with the two-gap Lame 
potential 

It was shown in |T^ by direct substitution that the equations ( |1.4| ) have the 
following periodical solutions: 

gi = Cdn(ax, k)sn{ax, k), q2 = Cdn{ax, k)cn{ax, k), (4.3) 

where sn, cn, dn are the standard Jacobian elliptic functions 0, /c is the 
modulus of the elliptic functions < k < 1, and the characteristic parameters 
of the wave: amplitude C, temporal pulsewidth 1/a and k are related to the 
physical parameters Q and, Qq through the following dispersion relations 



C2 = ^ + 21^0, k^ 



lOl^o 



a 



15 ■ 



(4.4) 



Another elliptic solution was found in |pO| . 

qi = Ciakcn{ax, k)sn{ax, k),q2 = C2sn^{ax, k) + C3, (4.5) 

where C, Ci, a and k are expressed through physical parameters VL and VLq 
by the following relations 



a 



6c 

2 



1 + 



cl 



15 



-'1 e 

10' 



20c 



(c- 15)(c + 5)' 



(4.6) 



where (? = 15Q/Qq 
The solutions 



I3| J4.5| ) can be obtained from the general periodic solution 
T^j. The solutions ([4. 3]) have the following spectral interpretation. 



They are the Lame polynomials associated with the eigenvalues 3e2, 3e3 (see 
Table 0), 

6 



ei 
6 



63 



-{p{x + Us) - ei){p{x + Us) - 63) 



-{p{x + UJ3) - ei){p{x + UJ3) 



62 



ei - 63 

One can show that the relations (p.ll|) become ( [4.4|) . Analogously one can 
prove that the normalised Lame polynomials associated with the eigenvalues 
3e2, —a/3 (72 gives the solution ( |4.5| ). We prolong this analogue by constructing 
some new periodic solutions corresponding to the Treibich-Verdier potentials 
given in Table |I]. 
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4.3 Periodic solutions associated with the two-gap Treibich- 
Verdier potentials 

Below we construct the two periodic solutions associated with the Treibich- 
Verdier potential W3 in the same way as the solutions (|4.3|) and ( [4 .51) are 
connected with the two-gap Lame potential W4. Let us consider the potential 

U,{x) = 6p{x + u,) + ai^LZ^^li^i^ (4.7) 

p{x + UJ3) - ei 

and construct the solution in terms of Lame polynomials associated with the 
eigenvalues 2:1,2:2, zi > Z2 



Zl = 62 + 2ei + 2^/(61 - 62) (7ei + 2e 



Z2 = e3 + 2ei + 2^(ei-e3)(7ei + 2e3). (4. 
The finite and real solutions qi , q2 have the form 

iC p{x + uj^) - 63 (^^J {p{x + LU^) - ei 



ei- 62 + \J (ei - e2)(7ei + 262) 
?>^{p{x + uj3) -ei 



(4.9) 



q2 = C\l p{x + uj'i)-e2[\l (p(x + u^) - d 



^ ei - 63 + ^(ei - e3)(7ei + 263) \ 
3\/(p(a; + W3) - ei / 

Here to'3 is the pure imaginary half period, < t < 2co' with 2uj being the real 
period, and 

= > 0. (4.11) 

Zl - Z2 

The parameters fl and are linked with the Weierstrass parameters Cj 
through 

2fi = -5(^1 + Z2), 2^0 = zi-Z2 (4.12) 

with Zl, Z2 given by ( [4.8| ). By eliminating Cj from these formula and from the 
formula for the amplitude ([4.11|) we arrive at the dispersion relations, 

c^rio = 9, o<k'^ <i, 
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i 60 ^ —— . (4.13) 



5^0 P + 2v/r^Fv/4^ - 2 V4 - 3P ' " 15' 

Analogously we can find the elliptic solution associated with the eigen- 
values 



We have 



Zl = e2 + 2ei + 2^(ei-e2)(7ei + 2e2), Z2 = -Qei, (4.14) 

qi = C(p(x + cu3) -ei), (4.15) 
q2 = iC^J p{x + UJ3) - 63 (^^J {p{x + UJ3) - ei 



ei - 62 + A/(ei - e2)(7ei + 262) , 

+ / I , (4.16) 

where C is given by ( [4.11| ) but with zi, Z2 given by ( [4.14| ). The corresponding 
dispersion relation has the form 

C^no = 9, - < < 1, 
8 

< — < -. (4.17) 



5^0 5 - 2A;2 + 2^l^^A^ " Sf^o ~ 3 

The Treibich-Verdier potential also yields elliptic solutions, but the so- 
lutions corresponding to the potential given in the Table blow up. From 
general considerations we conjecture that there exist non blow-up real solu- 
tions associated with the isospectral deformation of this potential. 

The solutions derived represent stationary periodical waves. However, 
taking into account the invariance of the equations ( |1 . 1| ) under a Galilean 
transformation |l^, P], they also represent travelling periodical waves. 



5 Exact quasi-periodic solutions 

In this section we give the theta functional expressions for the trajectories of 
the system under consideration using the Rosenhain ultraelliptic functions 
[ p^ (see also [0), i.e. Abelian functions associated with an algebraic curve of 
genus two. We also show how to reduce ultraelliptic solutions to the elliptic 
solutions discussed before. 
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5.1 Theta functional integration 

Let {za,Z/^) G (zi,...,Z5) be two arbitrary branching points of the curve 
( |2.5|) . After the transformation {w,z) {C^O^ 

W = 2i{Zf3 - Za)\l{Zi - Za){Zj - Za){Zk " ^a)C, 
Z = izp-Za)^ + Za (5.1) 



the curve becomes 28 



= e(i - 0(1 - '^'0(1 - A^O(i - /^'O (5.2) 

with 

2 Zi3 — Za .2 Zj3 — Za 2 Zfj — Z^ /r- o\ 

K = , A = — , /i = — . (5.3) 

Zi Zq^ Zj Zq, Zj^ Zfy 



Let us fix on (p.5|) the homology basis, (A; B) = (Ai, A2; Bi, B2) and the 



conjugated basis of differentials of the first kind 

/ N Ciii + C12 ,^ C2ii + C22 ,^ 
V = [vi, V2), vi = d4 V2 = 

normalised as 

(/ V, / v; / V, / v) = (l2;r), (5.4) 
\JAi Ja2 Jbi Jb2 ) 

where I2 is the 2x2 unit matrix and the period matrix r belongs to the 
Siegel upper half space ^2 of degree 2. The Jacobi inversion problem can be 
rewritten as 



■wi + / -ui = 2cii J^^j^-^x = MiX + Mio, (5.5) 



MO MO 
Ml /■M2 



/ V2+ [ V2 = 2c2iW = U2X + U20 (5.6) 

J fj,0 J fio V KAU 



The solution of the problem is expessed in terms of theta functions with 
characteristics [e] defined on x S2 



e 



'e' ' 




e" 


(w 



e' , e'A . . e' e" 



J2 exp jvri (^{n + -)r, {n + -)j + 2m {n + -,w + — 



(5.7) 
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where (, ) is a Euclidean scalar product and Im r is positive definite. For 
integer characteristics we have 



6 



'e' ■ 




e" 


(w 



expvrz 



X w 



\{e',Te') + {e'^) + \{e\e") 




(5.8) 



where 9{z\t) = 9 [q] {z\t). The characteristic [e\ is called even if the corre- 
sponding theta function is an even function {{e,e') = (mod2)) and odd if 
the corresponding theta function is an odd function {{e,e') = 1 (mod2)). 

The values of theta functions and their derivatives at the point z = are 
called theta constants which we denote by ^[e:](0;r) = 9[e], if the character- 
istic [e] is even and d9[e]{z] T)/dzi\^^Q = 9i[e\,i = 1,2 if the chacteristic [e] is 
odd. 

The function ( |5.7| ) satisfies two sets of functional equations , the trans- 
formational property 

^[£](w + n" + nV|r) 

= exp7rz[-(nV,n') -2(n",w) + (e', n') - (£",n')] ^[£](w|r) (5.9) 

where n', n" e and the modular property, which describes the transforma- 
tion of the theta function under the action of the group 5*^4 (Z). 

The branching points are expressed in terms of theta constants as 



0' m n ' 



0' [J?] [^^] 



9^ [g?] 



[^a [^g] 
0' [g?] 



The Jacobi inversion 


problem 


tias the solution 


0' [1?] 




> - 


Xq) + Uo 


r) 


= -kA/I yUl/i2, 


9^ [??] 


(u 




Xo) + Uo 


r) 




(u 


'x — 


Xo) + Uo 


r) 


/tA/i (1 - /ti)(l - /i2), 


9^ [??] 


(u 


[x- 


Xo) + Uo 


r) 


VI - kVi - AVI - /i' 


[SI] 


(u 


[x — 


Xo) + Uo 


r) 


A/i (1 — /t^/ii)(l — K^/i2) 


9^ [??] 


(u 


[x- 


Xo) + Uo 


r) 


Vl — — X^\J — /i^ 




(u 


[x — 


Xo) + Uo 


r) 


K/i (1 - A^/ii)(l - A^/i2) 


9^ [??] 


(u 


[x- 


Xo) + Uo 


r) 


VI - A VA' - /i V/t^ - A2 


n 


(u 


[x - 


Xo) + Uo 


r) 


kA (1 -/i2^i)(l -/i^/is) 


9- [??] 




[x- 


Xo) + Uo 


r) 


VI - /iV«^' - /iVA^ - /i^ 



(5.10) 
(5.11) 
(5.12) 
(5.13) 
(5.14) 
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The solutions of equations ( |1.4| ) are expressed in terms of the formulae 
(|5.10| - ^.14D . To do this we have to put -21,^2 into correspondence with two 
branching points of the curve ( |5.2| ). In particular, let us choose Zi^2 in such 
a way that their images Zi^2 under the transformation (|5.1| ) becomes the 
branching points and and the images of /ii,2 niove inside the gaps 
^2 = < /i2 < /i^, 5i = < /ii < 1. Let us also fix the vector uq as a half 
period uq = (l/2)(rii + T12, T12 + r22 + 1), which shifts the characteristics at 
Then we find the following quasiperiodic solutions 



2 _ , [Zp - Zj ^v^r^y/A^^y/;?^ 



[Zx - Z2)KjJ, 

[Ig] (u(x-xo) +Uo|r) 

^2 _ 2 (^^/^-^g)' 9^l\] {u{x - Xq) + Uo\t) 
^ {Zi - Z2)KXfJ^6^[ll]iu{x - Xq) +uo\t)' 

where the frequences Ui^2 are in general noncommensurable, 

^ - Za62 [qI] ^ - ^Q^l [qI] (^^7\ 

We emphasise that the solution given is parametrised only by the the r- 
matrix and expressed in terms of theta functions and theta constants which 
are rapidly convergent and therefore this expression is convenient for numer- 
ical calculations. 



5.2 Reduction to elliptic functions 

The given solutions are quasiperiodical functions of time. They are reduced 
to elliptic functions under some restrictions on the parameters of the system 
which can be formulated in terms of conditions on the matrix r. 

More precisely, it follows from the explicit formulae (|5.15| , |5?T6|) for qi , q2 
that the following two conditions are sufficient to require the solution be an 
elliptic function of x. 

1. The matrix r is taken in the form 

r=\^\' iVGN (5.18) 

\ N ^22 y 
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2. Ml or U2 = 0. 

It follows immediately from the transformational properties of theta func- 
tions that under these conditions the solution is an elliptic function of the 
A^-th order. The condition 2 is reduced with the help of the formulae ( |5.17| ) to 
the condition of the vanishing of some theta constants. It is remarkable that 
the conditions 1, 2 are also sufficient conditions, as was proved in [|T1| within 
the Weierstrass reduction theory of theta functions to lower genera 0, ^ . If 
only the condition 1 is satisfied then the solution (|5.15| , |5n6|) is quasiperiodic 
and expressed in term of two elliptic functions with the Jacobian moduli 
A^Tii and Nt22- 

To demonstrate how this approach works we derive the elliptic solution 
(|4.9| , |4.10|) from the theta functional approach. To this end we fix the period 
matrix in the form ( |5.18D and put = 4. We remark that this case of the 
reduction of theta functions was studied by Bolza (see also 0). 

The computation is based on the addition theorem of the second order 
(see e.g. 

(5.19) 



e[e]{^\T)e[S]{y\T) 



e" + 6" 



(x + y|2r)e 



-" - 5" 



e 



where the summation runs over p = (0, 0), (0, 1), (1, 0), 



(x 

1, 1) and 



y|2r) 



e 





(5.20) 



-1 ^2 
^1 ^ ^2 ^2^ ^1 

We now show that under conditions 1° and 2° the argument of the theta 
function in ( p.l0fj5.1^ ) becomes [xjluj^ 0). The condition M2 = or Q\ [J J] = 
becomes the condition (see Appendix) 



(5.21) 



where i^i = ^9.(0;4rii), = ^?i(0;4r22) 



2, 3, 4. The condition ( |OT[ ) is 



equivalent to the relations between Jacobi theta constants 



1-4 




-4 



(5.22) 
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Using the explicit expressions for the branching points (see Table 1) we find 
after substitution of the corresponding theta constants (see Appendix) in 

Ui = — ^ — , 



Chosing Z/3^a = 63 + 2ei ± 2^ (ei — 63) (7ei + 263) we find after the transfor- 
mation E = BAB (see e.g.p]), Ui = l/2uj = a/ci — e3/7i{}l, in accordance 
with §. 

The same arguments permit us to compute ( |5.15| , pTT6D in terms of elliptic 
functions. Let us consider the theta function entering into the expression for 
( |5.15| ). Applying the addition theorem ( [5 .191 ) and putting W2 = according 
to condition 2° and using the expressions for theta constants we find 



i'ol] ( 



2uj '^1' 



0[Vi] (^,0|r 



^2 (^) 



31^ 



X 



Xo 



2uj 



C2 



Bl] ( 



x-XQ nl 
2w '"I' 



X 



) ^1 (^) 



2u 



+ 



X — Xq 

2uo 



2uj 



+ 



2„q2 fx — Xq 



X 



Xq 



2uj 



X 



Xo 



2u 



(5.23) 



(5.24) 



where Ci,2 can be computed using the theta constants from the Appendix, 



Co 



5/2 



(5.25) 



where s = y^j — 4'j^4, t = yi^i —4^92 One can see that substituting ( |5.23| , 
5.24 ) and ( p. 15 , 3. 14 ) performing the transformation BAB (see, e.g. 0]) 
and setting xq = we obtain the elliptic solutions ( [4.9| , [4.10| ). In the same 
way the another elliptic solution associated with two-gap Treibich-Verdier 
potential ( |4.15|J4.1^ ) can be obtained from the theta functional solutions. 
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A Theta Constants 



We denote the Jacobi theta constants by ■dj — ■dj (0|2^tii), ■dj — ■dj (0| 2^X22), 
J = 2, 3, 4. 



p=l: Let T 



Til 

1 

2 

10' 
01 



1 

2 

T22 



Then 



(2^?2^3^3^4)'/', e 



01 
01 
01 

01 
01 



= -770 

— —m0 

= —1710 

— 17:0 

= -7:0 

= -7:0 



02 



11 

00 




02 


11 
01 


"01' 

00 


'^2! 


02 


"01" 
01 



01' 


^2) 




"01" 


10 


02 


11 



'10' 
01 


1/4, 


02 


"10" 
11 


'lO' 
00 


V4, 


02 


"10" 
10 



= —17:0 
= -7:9 
= -7:9 
-- -7:9 
— 17:0 
= —ir:0 



1/4, 



1/4, 



7?^ 



p=2: Let r = 



^ 1 

ni 4 

1 ^ 

I T-22 



and denote X = ^^■d^, Y = 'i?2i?4, ^ = ^i'&2, 



A = -X^ + Y^ + Z^, B = X^-Y^ + Z"^, C = X^ + Y^-Z^, D = A + B + C. 
Then the following formulae hold: 




0' 
0' 
0' 



X + Y + Z, 



X + Y-Z, 



= X -Y + Z, 9 [??] =X -Y-Z, 

= 2V2(xy)V2(DV2 + ^z), 0^ ' 

= 2V2(XZ)V2(dV2 + ^Y), 9^ 
= 2V2(rZ)i/2pi/2 ^ ^^^^ ^2 



: 2V2(xy)V2 (2)1/2 _ 

: 2^/2(XZ)V2(dV2 _ ^y)^ 

2V2(rZ)l/2pl/2 _ 
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"10" 
10 




62 


10' 
10 


= -in{2XYy/\^lB'/^ + V2^lZ){D^/^ + V2Z)-'/\ 




01' 

01 






"01" 

01 


= -n{2X Zy/\^lC'/' + V2^lY){D'/^ + V2Y)-'/\ 




11' 

01 


= -t7ri2ZYy/\^lC'/^ + V2^lX)iD'/^ + v^X)-^^ 


02 


11" 

01 


= -7r(2Zr) V4(^2^i/2 + v^^2^)(L»V2 + V2X)-V2, 


01 


ll' 

10 


= -7r(2Zy)V4(^2^V2 + ^^lx){D'/^ + V2X)-^/^, 


02 


11" 

10 


= -7ri2ZYy/\^lB^/^ + V2^lX){D'/^ + V2X)-'/^, 


01 


"10" 
11 


= -i7i{2XYy/\^lA'/^ - V2z^lZ){D'/^ + V2Z)-'/^ 


02 


10' 
11 


= -i7T{2XY)'/\^lA'/^ - V2^lZ){D^'^ + V2Z)-^/^ 


01 


01' 
11 


= 7r(2XZ)V4(^2^V2 _ ^i^2y^(^l/2 ^ ^y)-i/2^ 


02 


01" 
11 


= -i7r(2XZ)l/^(l?2AV2 _ v^i^2y)pl/2 ^ ^y)-i/2 



References 

[1] H Airault, H P McKean, and J Moser. Rational and elliptic solutions of 
the KdV equation and a related many-body problem. Comm. Pure and 
Appl. Math., 30:94-148, 1977. 

[2] H Bateman and A Erdelyi. Higher Transcendental Functions. Volume 
2, McGraw-Hill, New York, 1955. 

[3] E D Belokolos, A I Bobenko, V Z Enolskii, A R Its, and V B Matveev. 
Algebraic-geometrical Methods in the Theory of Integrable Equations. 
Springer, Berlin, 1994. 

[4] E D Belokolos and V Z Enolskii. Isospectral deformations of elliptic 
potentials. Russian Math. Surveys, 44:5:155-156, 1989. 

[5] K J Blow, N J Doran, and D Wood. Polarization instabilities for solitons 
in birefringent fibers. Opt. Lett, 12:202-204, 1987. 



21 



O Bolza. Ueber die Reduction hyperelliptischer Integrale erster Ordnung 
und erster Gattung auf elliptische durch eine Transformation vierten 
Grades. Math. Ann., 28:47-496, 1887. 

D N Christodoulides. Black and white vector solitons in weakly bire- 
fringent optical fibers. Phys.Lett.A, 132:451-452, 1988. 

D N Christodoulides and R I Joseph. Vector solitons in birefringent 
nonlinear dispersive media. Opt. Lett., 13:53-55, 1988. 

R J Dowling. Stability of solitary waves in a nonlinear birefringent 
optical fiber. Physical Review A, 42:5553-5560, 1990. 

J C Eilbeck and V Z Enolskii. Elliptic Baker-Akhiezer functions and 
an application to an integrable dynamical system. J. Math. Phys., 
35(3):1192-1201, 1994. 

J C Eilbeck, V Z Enolskii, V B Kuznetsov, and D V Leykin. Lin- 
ear r-matrix algebra for systems separable in parabolic coordinates. 
Phys.Lett.A, 180:208-214, 1993. 

J C Eilbeck, V Z Enolskii, V B Kuznetsov, and A V Tsiganov. Linear 
r-matrix algebra for classical separable systems. J. Phys. A, 27:567-578, 
1994. 

J C Eilbeck, P S Lomdahl, and A C Scott. The discrete self-trapping 
equation. Physica D, 16:318-338, 1985. 

V Z Enolskii and J C Eilbeck. On the two-gap locus for the elliptic 
Calogero-Moser model. 1994. Preprint |hep-th/9404042| . 

V Z Enolskii and N A Kostov. On the geometry of elliptic solitons. Acta 
AppUcandae Matheniaticae, 1994. 

M Florjanczyk and R Tremblay. Periodic and solitary waves in bimodal 
optical fibers. Phys.Lett.A, 141:34-36, 1989. 

J Hietarinta. Direct method for the search of the second invariant. 
Physics Reports, 147:87-154, 1987. 



22 



[18] M F Jorgensen and P L Christiansen. Hamiltonian structure for a mod- 
ified Discrete Self-trapping Equation. Chaos, Solitons, and Fractals, 
4:217-225, 1994. 

[19] N A Kostov. Quasi-periodical solutions of the integrable dynamical 
systems related to Hill's equation. Lett.Math.Phys., 17:95-104, 1989. 

[20] N A Kostov and M Uzunov. New kinds of periodical waves in birefringent 
optical fibers. Opt.Commun., 89:389-392, 1992. 

[21] A Krazer and W Wirtinger. Abelsche Funktionen und allgemeine Theta- 
functionen, pages 603-882. Teiibner, 1915. 

[22] N Manganaro and D F Parker. Similarity reductions for variable- 
coefficient coupled nonlinear Schrodinger-equations. J. Phys. A, 26: 
4093-4106, 1993. 

[23] H P McKean and P van Mocrbeke. The spectrum of Hill's equation. 
Invent. Math, 30:217-274, 1977. 

[24] C R Menyuk. Nonhnear pulse-propagation in birefringent optical fibers. 
IEEE J.Quan.Electron., 23:174-176, 1988. 

[25] D Mumford. Tata lectures on theta, Vol.1, Vol.2. Birkhauser, Boston, 
1983, 1984. 

[26] S P Novikov. The periodic problem for Korteweg de Vries equation. 
Funct. Analiz i ego Prilozen., 8:54-66, 1974. 

[27] V Ravoson, A Ramani, and B Gramaticos. Generalized separability for 
a Hamiltonian with nonscparablc quartic potential. Phys. Lett. A, 1994. 

[28] G Rosenhain. Abhandlung uber die Funktionen zweier Variabler mit 
vier Perioden. Mem. pres. I'Acad de Sci. de France, des savants, IX:361- 
455, 1851. 

[29] R Sahadevan, K M Tamizhmani, and M Lakshmanan. Painleve analysis 
and integrability of coupled non-linear Schrodinger-equations. J. Phys. A, 
19:1783-1791, 1986. 



23 



[30] M V Tratnik and J E Sipe. Bound solitary waves in a birefringent optical 
fiber. Phys.Rev.A, 38:2011-2017, 1988. 



[31] A Treibich and J-L Verdicr. Solitons elliptiqucs. In Special volume for 
60th. anniver. of Prof . A.Grothendieck, Birkhaiiser, Boston, 1991. 

[32] V E Zakharov, S V Manakov, S P Novikov, and L P Pitaevskii. Soliton 
theory: inverse scattering method. Nauka, Moscow, 1980. 

[33] V E Zakharov and E I Schulman. On the integrabihty of the system of 2 
coupled non-hnear Schrodinger-equations. Physica D, 4:270-274, 1982. 



24 



